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00 ! Abstract 

For several processes at colliding beams, macroscopically large impact param- 
, eters give an essential contribution to the standard cross section. These impact 

\^ ' parameters may be much larger than the transverse sizes of the colliding bunches. 

■ In that case, the standard calculations have to be essentially modify. The corre- 

, sponding formulae for such a beam-size effect were given twenty years ago without 

I taking into account correlations of particle coordinates in the beams. In the present 

' paper we derive formulae which necessary to take into account quantitatively the 

effect of particle correlations in the spectrum of bremsstrahlung as well as in pair 
\ production. Besides, we consider critically recent papers |17| llSj in which it was 

calculated a new additional "subtraction term" related to the coherent contribution 



Ph. into beam-size effect. We show that this result is groundless and point out the origin 

of the mistake. 



i INTRODUCTION 

The so called beam-size or MD-effect is a phenomenon discovered in experiments ^ at 
the MD-1 detector (the VEPP-4 accelerator with e^e~ colliding beams , Novosibirsk 
1981). It was found out that for ordinary bremsstrahlung, macroscopically large impact 
parameters should be taken into consideration. These impact parameters may be much 
larger than the transverse sizes of the interacting particle bunches. In that case, the 
standard calculations, which do not take into account this fact, will give incorrect results. 
The detailed description of the MD-effect can be found in review 0. 

In 1980-1981 a dedicated study of the process e^e~ — > e"'"e~7 has been performed 
at the collider VEPP-4 in Novosibirsk using the detector MD-1 for an energy of the 
electron and positron beams = Ep = 1.8 GeV and in a wide interval of the photon 
energy E^ from 0.5 MeV to E^ k, E^.. It was observed P that the number of measured 
photons was smaller than that expected. The deviation from the standard calculation 
reached 30% in the region of small photon energies and vanished for large energies of the 
photons. Yu.A. Tikhonov j3] pointed out that those impact parameters which give an 
essential contribution to the standard cross section, reach values of ~ 5 cm whereas 
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the transverse size of the bunch is a_[_ ~ 10~^ cm. The hmitation of the impact parameters 
to values ^ < cr^ is just the reason for the decreasing number of observed photons. 

The first calculations of this effect have been performed in Refs. ^ and [3] using 
different versions of quasi-classical calculations in the region of large impact parameters. 
Further experiments, including the measurement of the radiation probability as function 
of the beam parameters, supported the concept that the effect arises from the limitation 
of the impact parameters. Later on, the effect of limited impact parameters was taken 
into account when the single bremsstrahlung was used for measuring the luminosity at 
the VEPP-4 collider and at the LEP-I collider 0. In the case of the VEPP-4 experi- 
ment [HI , it was checked that the luminosities, obtained using either this process or using 
other reactions (such as the double bremsstrahlung process e^e~ —>■ e"'"e~77, where the 
MD-effect is absent), agreed with each other. 

A general scheme to calculate the finite beam size effect had been developed in pa- 
per jH] starting from the quantum description of collisions as an interaction of wave pack- 
ets forming bunches. Since the effect under discussion is dominated by small momentum 
transfer, the general formulae can be considerably simplified. The corresponding approx- 
imate formulae were given in [S]. In a second step, the transverse motion of the particles 
in the beams can be neglected. The less exact (but simpler) formulae, which are then 
found, correspond to the results of Refs. j3] and It has also been shown that similar 
effects have to be expected for several other reactions such as bremsstrahlung for colliding 
ep-beams jSj, [Ej, e+e"- pair production in e^e and •ye collisions [Hj. The correspond- 
ing corrections to the standard formulae are now included in programs for simulation of 
events at linear colliders. The influence of MD-effect on polarization had been considered 
in Ref. |TT]. In 1995 the MD-effect was experimentally observed at the electron-proton 
collider HERA J2] on the level predicted in [TUj . 

The possibility to create high-energy colliding beams is now wildly discussed. 

For several processes at such colliders a new type of beam-size effect will take place — the 
so called linear beam-size effect [13j. The calculation of this effect had been performed 
by method developed for MD-effect in [S]. 

It was realized in last years that MD-effect in bremsstrahlung plays important role for 
the problem of beam lifetime. At storage rings TRISTAN and LEP-I, the process of a 
single bremsstrahlung was the dominant mechanism for the particle losses in beams. If 
electron loses more than 1 % of its energy, it leaves the beam. Since MD-effect reduced 
considerable the effective cross section of this process, the calculated beam lifetime in 
these storage rings was larger by about 25 % for TRISTAN 14J and by about 40 % for 
LEP-I ^5] (in accordance with the experimental data) then without taken into account 
the MD-effect. According to our calculations at B-factories PEP-II and KEKB the 
MD-effect reduces beam losses due to bremsstrahlung by about 20%. 

It is seen from this brief listing that the MD-effect is a phenomenon interesting from 
the theoretical point of view and important from the experimental point of view. In the 
present paper we consider once again the MD-effect having in mind two aims. The main 
purpose is to take into account correlations of particle coordinates in the beams. Usually 
these correlations are small, however, more accurate measurement may be sensitive to 
them. In the present paper we derive formulae which necessary to take into account 
quantitatively the effect of particle correlations in the spectrum of bremsstrahlung as well 
as in pair production. 

Besides, we would like to consider critically recent papers [T71 HH] in which previous 
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Figure 1: Block diagram of radiation by the electron. 



results in 13 El about bremsstrahlung spectrum had been revised. It was claimed that 
an additional subtraction related to the coherent contribution has to be done. Analy- 
sis, performed in paper jT2j, results in the conclusion that this additional "subtraction 
term" in the spectrum is not important for the MD-1 experiment jl], but it should be 
taken into account in processing the HERA experiment ^2]; it also may be important 
for the future experiments at linear e+e" colhders. It should be noted that in paper [T7j 
there is no derivation of the starting formulae. On the other hand, in paper ^7] there 
is a general remark that their consideration was motivated by corresponding calculations 
for bremsstrahlung of ultra-relativistic electrons on oriented crystals. In our critical re- 
mark |16j we had shown that the starting formulae of J7] are incorrect. After that new 
paper [TH] was appeared in which there is the "derivation" of the starting formulae used 
in Unfortunately, this derivation is incorrect as well. 

In the present paper we analyze the coherent and incoherent contributions in the 
conditions, considered in papers ^3 ^] , when the coherent length /coh is much smaller 
than the bunch length / but much larger than the mean distance between particles a, i.e. 
o hoh ^ We derive expressions for the coherent and incoherent contributions and 
show that under these conditions the coherent contribution is completely negligible and, 
therefore, there is no need to revise the previous results. This conclusion is quite natural. 
A usual bunch at colliders can be considered as a gaseous media with a smooth particle 
distribution which has characteristic scales of the order of bunch sizes. In particular, the 
average particle density in such a bunch has the only scale in the longitudinal direction 
— the length of the bunch /. Therefore, the average field of the bunch has the spectral 
components in the region of frequencies uj = q^c ~ c//coh ^ c// and vanishes in the region 
of much higher frequencies considered here. On the contrary, in the crystal case there is 
another scale related to the size of the particle localization in the crystal structure. In this 
case, the additional subtraction should be taken into account for incoherent contribution. 
It seems that the electron radiation on oriented crystals played a misleading role for 
consideration of the MD-effect in [T7|IIH1- To clarify a question we give our calculations 
in full details and pointed out the origin of the corresponding mistake in [T71 ITH]. 

In next Section we present the qualitative description of the MD-effect. In Sect. Ill 
we discuss our approximations. Basic formulae for coherent and incoherent contributions 
are given in Sect. IV. Correction to the standard cross section, related to the particle 
correlations, are derived in Sect. V. Our critical remarks about results of papers [T7| ITSj 
are presented in Sect. VI. Some conclusions are given in Sect. VII. 
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Figure 2: Compton scattering of equivalent photon on the electron. 

2 QUALITATIVE DESCRIPTION OF THE MD-EFFECT 

Qualitatively we describe the MD-effect using as an example the ep — *■ epj process^. This 
reaction is defined by the diagrams of Fig. ^ which describe the radiation of the photon 
by the electron (the contribution of the photon radiation by the proton can be neglected). 
The large impact parameters Q ^ cr^, where a± is the transverse beam size, correspond 
to small momentum transfer hq± ~ {^/q) ^ In this region, the given reaction 

can be represented as a Compton scattering (Fig. ^ of the equivalent photon, radiated 
by the proton, on the electron. The equivalent photons with frequency u form a "disk" 
of radius Qm ~ 'Jpc/u where = Ep/{mpC^) is the Lorentz-factor of the proton. Indeed, 
the electromagnetic field of the proton is 7p-times contracted in the direction of motion. 
Therefore, at distance q from the axis of motion a characteristic longitudinal length of a 
region occupied by the field can be estimated as A ~ ^/7p which leads to the frequency 
~ c/A ~ 7pc/ Q. 

In the reference frame connected with the collider, the equivalent photon with energy 
tijjj and the electron with energy E^. ^ huj move toward each other (Fig. E)) and perform 
a Compton scattering. The characteristics of this process are well known. The main 
contribution to the Compton scattering is given by the region where the scattered photons 
fiy in a direction opposite to that of the initial photons. For such a backward scattering, 
the energy of the equivalent photon hw and the energy of the final photon and its 
emission angle 9^ are related by 

and, therefore, for typical emission angles ^7 ^ l/7e one has 

As a result, we find the radius of the "disk" of equivalent photons with the frequency 
u (corresponding to a final photon with energy E^) as follows: 

f^n. = ^ ~ A, 47e7p ^^V^ ' A, = -^ = 3.86 -10-^1 cm. (3) 

uj E^ rrieC 



^Below we use the following notations: N^. and Np are the numbers of electrons and protons (positrons) 
in the bunches, az — I is the longitudinal, (7^ and ay are the horizontal and vertical transverse sizes of the 
proton (positron) bunch, 7^ = E^l {mf-c^), 7p — Ep/^nipC^) and ^ e'^ / (rrieC^) is the classical electron 
radius. 
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For the HERA collider with Ep = 820 GeV and = 28 GeV one obtains 

gm>lcm for < 0.2 GeV . (4) 

Equation (jH)) is also valid for the e~e~^ e~e"'"7 process with replacement protons by 
positrons. For the VEPP-4 collider it leads to 

^„ > 1 cm for E^ < 15 MeV , (5) 

for the PEP-II and KEKB colliders we have 



> 1 cm for E^ < 0.1 GeV . (6) 

The standard calculation corresponds to the interaction of the photons forming the 
"disk" with the unbounded flux of electrons. However, the particle beams at the HERA 
collider have finite transverse beam sizes of the order of a± ~ 10~^ cm. Therefore, the 
equivalent photons from the region cr± ^ g ^ Qm cannot interact with the electrons from 
the other beam. This leads to the decreasing number of the observed photons and the 
"observed cross section" daohs is smaller than the standard cross section da calculated for 
an infinite transverse extension of the electron beam, 

daohs = da - da^ov (7) 
Here the correction dacov can be presented in the form 

dacor = daciyj, E^, E^) dn{u) (8) 

where dn{uj) denotes the number of "missing" equivalent photons and dac is the cross 
section of the Compton scattering. Let us stress that the equivalent photon approximation 
in this region has a high accuracy (the neglected terms are of the order of l/7p). But for 
the qualitative description it is sufficient to use the logarithmic approximation in which 
this number is (see ^20j, §99) 

a duj dq"^, , , 

dn = ^ . 9 

71 to qj_ 

Since q± ^ 1/ g, we can present the number of "missing" equivalent photons in the form 

dn=^^^ (10) 
7C to g'^ 




Figure 3: Scattering of equivalent photons, forming the "disk" with radius g^ = Ipc/to, 
on the electron beam with radius a±. 
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with the integration region in q: 



a±<Q<Qm = ^. (11) 



As a result, this number is equal to 



dnM = 2^^1n^, (12) 



and the correction to the standard cross section with logarithmic accuracy is 

rf^c, = ^ (1 - !/ + m in ^J£!£^^ (13) 

3 y ^ ^ ' y(J±_ Ee 

3 APPROXIMATIONS 

For future linear e+e~ colliders the transverse sizes of the beams will change significantly 
during the time of interaction due to a mutual attraction of very dense beams. However, 
for most of the ordinary accelerators, including practically all e+e~ and ep storage rings, 
the change of the transverse beam sizes during the collisions can be neglected. Below 
we use two main approximations: 1) the particle movement in the bunches has a quasi- 
classical character; 2) the particle distribution remains practically unchanged during the 
collision. For definiteness, we use again the ep collisions as an example. 

Therefore, if the proton (electron) bunch moves along (opposite) the direction of 2;-axis 
with the velocity Vp (vg), its density has the form 



n. 



^ Upig, Z - Vpt) , Ue ^ ne{Q, Z + Vet) . (14) 



We also introduce so called "transverse densities" 



npio) = J Updz, Ueio) = J Uedz (15) 

which is equal to the total number of protons (electrons) which cross a unit area around 
the impact parameter g during the collision. Using the transverse densities, we express 
the luminosity for collisions of beams whose axes separated by impact parameter g as 

Hq) = J ne(r±) np{g + r_L) rf^r^ . (16) 

The usual luminosity for a single collision of ep-beams Lgp is then 

Lep = m . (17) 



^Within this approximation, the standard cross section has the form 

da = dac ^ = —ari — [1 - y + |r) In 

with the integration region Lo/{cjp) < q'x < rriec/h corresponding to the impact parameters q in the 
interval Ag ^ g ^ Qm- 
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Below we consider in detail the case when an electron deflection angle Of, is smaller 
than the typical radiation angle ~ l/7e. It is easy to estimate the ratio of these angles. 
The electric E and magnetic B fields of the proton bunch are approximately equal in 
magnitude, |E| ^ |B| ~ eNp/{ax + CTy). These fields are transverse and they deflect the 
electron into the same direction. In such fields the electron moves around a circumference 
of radius R ~ 7emeC^/(e-B) and gets the deflection angle Oe ~ IjR- Therefore, the ratio 
of these angles is of the order of 

6 r N 

77^-^ = — (18) 
(l/7e) a^ + ay 

The parameter ?7 ^ 1 only for the SLC and future linear e~^e~ colliders, in most of the 
colliders f] ^ 1- 

In our consideration we use the equivalent photon approximation. In the region of 
interest (where impact parameters are large, Q ^ o'±) this simple and transparent method 
has a high accuracy: the neglected terms are of the order of I/7. It should be stressed 
that the operator quasi-classical method, used in Ref. [T7j, just coincides in this region 
with the equivalent photon approximation. 



4 COHERENT AND INCOHERENT CONTRIBUTIONS 



4.1 General formulae 

The corresponding formulae for the number of events in a single collision of the electron 
and proton bunches can be found in papers ^IT], f^. To calculate the MD-effect, we need 
to know the distribution of equivalent photons (EP) for large values of impact parameters. 
In this region we can consider the electron-proton scattering as the scattering of electrons 
on the electromagnetic field of the proton bunch. Replacing this field by the flux of EP 
with some frequency distribution, we obtain the number of events in the form 

dN = dL^ei^^) daciuj , Ee, E^) , dL^^i^) = n^(Q,uj)dujne{Q) d'^g . (19) 

Here ^e(^) is the transverse electron density (fT3j) and n.y{Q,uj)duj is the transverse density 
of EP with the frequencies in the interval from u to uj + du. The quantity dL^^i^^) denotes 
the differential luminosity for the collisions of EP and electrons and dac{uj, Eg, E^) is the 
Compton cross section for the scattering of the equivalent photon with the frequency 00 
on the electron. 

For comparison with the experimental data the number of events in a single collision 
of beams dN should be averaged over many collisions of bunches in a given experiment. 
For example, the typical rate at the HERA collider is less then 1/100 bremsstrahlung 
photons in a certain interval of frequencies per a single collision of the beams, therefore, 
in that experiment the averaging over many collisions of bunches really does exist. 

The transverse density of the EP is determined by density of the electromagnetic field 
for a given frequency, i.e. by | ^uj{q) P /(4vr), where E^(^) is the spectral component of 
the collective electric field of the proton bunch. As a result, the transverse density of the 
EP is 

n^{Q,uj) duj=-^ {\ E^q) ^, (20) 
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where the sign (. . .) denotes the above mentioned statistical averaging. The field Euj{q) 
itself depends on a distribution of charges in the proton bunch at t = 0. We introduce 
the exact (fiuctuating) density of the proton bunch n{r) and the averaging density 

n,(r) = (n(r)) (21) 

as well as the corresponding form factor 

Fp(q) = j ripir) e~'^' d?r (22) 

with the normalization 

Fp(0) = j np{r) dh = Np . (23) 

In the classical limit 

n(r) = J]5(r-r,), (24) 

a 

where is the radius- vector of the a-th proton. In these notations, the exact (fiuctuating) 
collective field is 



(26) 



with = u/c. 

As a result, the number of events 

depends on the beam structure factor 

5(q, q') = j Sir, r') e^'^'''-''''''^ d\ d\' , ^(r, r') = ( n(r) n(r') ) (27) 

in which 

q, = q'^ = u/c. (28) 

Below we analyze these formulae in conditions when the coherence length /coh ^ c/uj 
is much smaller than the bunch length /, but much larger that the mean distance between 
particles in the beam a, i. e. at 

a « - ~ /coh = -^r- (1 - E,/E,) « /. (29) 
4.2 The beam structure factor 

The obtained general formulae include the coherent and incoherent contributions. The 
coherent contribution is determined by the average field which is given by Eq. ()25|) with 
the replacement of the exact density n(r) by the average density np(r). The averaged 
density of the proton bunch has a single scale in the longitudinal direction — the length 
of the bunch /. Therefore, the average field of the bunch is essential in the region of 
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frequencies u = cq^ ^ c/l and should be small in the region of large frequencies u 3> c/l. 
In particular, if the proton bunch has the Gaussian distribution, 



^p(r) = ,^^^.,/9. . , exp 



z'^ 



2al 2al 2P 



(30) 



(27r)3/V,a,/ 

its form factor is equal to 

Fp(q) = Np exp - i(gj,aj,)2 - ^{ul/cY] (31) 

and vanishes in the discussed region of frequencies from the interval ()29|). 
If we introduce the density fluctuation 

An(r) = n{r) - np{r) , (32) 

we can rewrite the average product of densities in the form 

{n{r)n{r')) = np{r)np{r') + (An(r) Ari(r')) . (33) 

In accordance with this presentation, we split the function ^(r, r') in two items called 
coherent and incoherent contribution: 

S = Scoh + 5'incoh , 5'coh(r, r') = np{r)np{r') , 5'incoh(r, r') = (An(r) An(r')) . (34) 
The coherent contribution to the structure factor is equal to 

5coh(q,q) = ^p(q)i^;(q)- (35) 



This formula was used in Refs. [21], [22] to obtain main characteristics of the coherent 
bremsstrahlung. It also allows us to obtain the following estimate for the Gaussian beam 
in the region of interest (at \qx\, \q'^\ ^ 1/cTx and \qy\, |g^| < l/o-y): 

^coh(q, q') ~ exp [-{col/cy] . (36) 

Let us now consider the incoherent contribution. A bunch at colliders can be treated as 
a continues media with a smooth average particle distribution. It was shown in Appendix 
that for such a media the function S'incoh(r, r') is expressed only via the average density 
rip(r) and via the correlation function C(r,r') as follows 

5incoh(r, r') = 5{r - r') ^^(r) + C(r, r') . (37) 

If we neglect the correlations of the particle coordinates in the beam, the correlation 
function C(r, r') vanishes, and we obtain (taking into account Eq. 



5'incoh(q,q') =^(q±-q'±)- (38) 

It is important that this expression is determined only by the transverse average density 
of the proton bunch and it does not depend on u. Formula (|HHj) has been used to derive 
the previous results about MD-effect (for details see review and Sect. V). For the 
Gaussian beam in the region of interest, we get from a useful estimate 



5incoh(q,q')~Arp. (39) 
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The correlations of the particle coordinates may arise due to Coulomb interaction of 
particles in the beam. In this case the characteristic quantity — the correlation length 
^corrci — is related to the Debay radius. It is evident that the correlations are negligible 
if the correlation length is much larger then the coherence length, i.e. at /correi ^ ^coh- 
According to an estimate [Hj it is just the case for the VEPP-4 experiment p. In any 
case, the correlations may give an essential correction to the standard bremsstrahlung 
cross section only if 

^correl ^ ^coh • 

(40) 

Therefore, the important quantity is the spectral component of the correlation function: 



C^(r±, rl) = j C(r, v') e-''^^'-''^/' dz dz' . 



(41) 



With this notation the incoherent contribution is now: 

^incoh(q, q') = i^(q± - ql) + ^(q^, ql) , (42) 

where 

CM±, q'±) = J CM±, r'^) e-'^^-^--^'^^'^^ d'r^ d'r'^ . (43) 

5 CORRECTION TO THE STANDARD BREMSSTRAHLUNG CROSS SEC- 
TION 

Let us compare the coherent and incoherent contributions for the Gaussian beams. In 
this case, the ratio 

^V^-^^^^^-^'^^'pH-'/c)] (44) 

is determined by the parameters ul/c. Since hw ~ E^/ [47^(1 - E^/E^)], it is also 
useful to introduced the coherence length ()29|) and the critical energy for the coherent 
bremsstrahlung 

£. = ^f^. (45) 

If the coherence length is large, /coh ^ /, or if the final photon energy is small, E^ < E^., 
the parameter ul/c ^ 1 and the coherent contribution is dominant. 

On the contrary, in the region of large photon energy, E^^ ^ E^, or small coherence 
length, /coh "C /, considered here, the incoherent contribution dominates. In particular, 
for Np ~ 10^^ the ratio dN''"^ / dN'"^^°^ is small even for ul/c = 6, 

jyycoh 

A^pe-3^«1, (46) 



and the coherent contribution becomes completely negligible. In this case the number of 
events for bremsstrahlung can be presented in the form (cf. (jZj)) 

^^incoh ^ ^^^^^ ^ ^^^^^ = da- da^ov , (47) 
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where Lgp is the luminosity (fT7|) of the ep-colhsions, da is the standard cross section for 
the ep — > ep7 process and dacor is the correction related to the MD-effect. Then we 
perform integration over and q'j^ using the well known equality 

^ d\^ = — ^ K^{g/b) (48) 



qi + (1/6)2 5 ^ 

where -ft'„(a;) denotes the modified Bessel function of third kind with integer index n 
(McDonald function). As a resul, we obtain the correction to the standard cross section 
in the form (cf. with the approximate formulae (jH)), (jl2j) ) 

dacor = daciuj, E^, E^) -— G{uj) , (49) 

where dac is the Compton cross section and the function G{uj) consists of two items 

G{uj) = G^^\uj) + G^^\uj) . (50) 

The first item represents the previous result for the MD-effect (without taking into account 
correlations) , 



r d^Q 




1 T^eL 


L ^(o)J 



Kl^QlQm), Qm = —, (51) 
UJ 



where L{q) is defined in ()16p. Some other useful expressions for G'^^^uj) as well as its 
asymptotics can be found in [2*. The second item is directly related to the correlation 
function P7I). (HT|l : 

G^'\uj) = _ f l^'lM c^(r^ + ri + e) K,{tJq^)K,{t'Jq^) d'r^d'r'^ . 

(52) 

Note that the main contribution to G^^\uj) (f^T|) is given by the region of large impact 
parameters (lllj) while the main contribution to G^'^\uj) is given by the region g ~ a±. 

The quantity dacduj/uj in (|49|) can be expressed via the energy E^ and the emission 
angle 9^ of the final photon as follows (taking into account relation (^) 

aduj ndy dz ^, . ^, . , , 1 + n , , 

dac = 2arl^——^^F{y,z), F{y, z) = 2{l - y)——-^ + y\ 53 

Ti UJ y (1 + zy [1 + zY 



where 



— , 2; = (^^7e) , re = — , 

Ee mgC 



y = -^, z = {e,j,y, r, = -. (54) 



6 CRITICAL REMARKS ABOUT RESULTS of Refs. |n|, [H] 

We derive the final expression for the incoherent contribution from general equations 
(fTIJ|) . and as a simple consequence of natural assumptions about the particle 
distribution in a proton bunch. It is useful to rewrite these equations in the form conve- 
nient for comparison with the corresponding equations in ^1 ^]. To do this, we note 
that the Compton cross section dac oc |eMcomptonP where eMcompton is the amplitude 
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of the Compton scattering for the EP with the polarization vector e = 'Ei^^[q)/ \ E(^(^) |. 
Therefore, the number of events in a given colhsions of beams is proportional to |Mp, 
where 

M = E^(e)Mcompton (55) 

is related to the probability amplitude of the process. Further, we use Eqs. fl^ and (jlHj) 
and present the collective field of the proton bunch Ei^(^) as a sum of fields of all protons: 

E^(e) = ^E(")(e), E(r)(^) = ^^iri(f>l/^Je— (56) 

where = ^ — is the impact parameter between the electron and the a-th proton and 
the parameter Qm = ^pcjuj is the radius of the "disc" of EP (see Fig. 3). 
As a consequence, the amplitude M is the sum 

M = J] m„ ^-'^'-^l^ , = — iTi {djQm) -^Mcompton , (57) 

where the item exp(—iu; 2:^/0) is the contribution to M related to the interaction of the 
electron with the a-th proton, while |Mp can be presented as the double sum 

|M|2 = ^ m„m* e-'^^"''-^'')/" . (58) 

We split this sum into sum with a = b and sum with a ^ b: 

|M|2 = Si + S2, Si = 5^|m,p, S2 = 5^m,m,*e-^-(^"-^'')/^ (59) 

Equations (j37j) — can be considered as the same starting formulae in our approach, 
based on the equivalent photon approximation, and in approach of f7| ^] , based on the 
operator quasi-classical method. However, further calculations are quite different. For 
simplicity, below we consider the case when we can neglect the correlations between the 
particle coordinates in the proton bunch. 

In our approach, the number of events is proportional to |Mp averaged over collisions 
of beams, i.e. 

rfAT oc (|M|2) = (Si) + (S2) (60) 

In the considered region of large frequencies (j29j) . the item (S2), corresponding to the 
coherent contribution, vanishes, 

1(^2)1 « (Si). (61) 

Since nia does not depend on the longitudinal coordinate Za, the average value of |map is 
determined by the transverse average density of the proton bunch fll5p . 

{\ma?) = J \ma\' ^ rfV, = J |m,p Sg^ , (62) 

and it does not depend on index a. Therefore, the item 

(Si) =iVp(|m,|2) (63) 
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leads to the correction, corresponding to C^^^ in (jlHl) — 

Authors of [T71 ITHI as the first step had performed averaging over transverse coordinates 
of the protons. Certainly, after that they get the same expression for Si as in ()63|) . For 
E2 they had obtained the following expression 



■\U}{Za-Zb)/c 



where 



ma 



N, 



d^Qa 



(64) 



(65) 



The principal mistake in ^Hj consists in the incorrect calculation of Z. It is not diffi- 
cult to understand its true behavior in the considered region of large frequencies (j^ . 
The quantity Z fluctuates near zero for various sets of coordinates {zg} = -21,^2, ■■■,ZNp, 
corresponding to various collisions of beams, and after averaging over many collisions 
one obtains the estimate (j61|) . This natural behavior of Z is confirmed by numerical 
calculations given below. 

When calculating Z, the authors of [TH] add and subtract the items with a = 6, as a 
consequence. 



J-Np, J 



a,b 



-iuj{Za-Zi,)/c 



-iujZa/c 



(66) 



Their next step consists in replacement the sum J by the integral 

2 



J 



-iujz/c 



np(r) (fr 



(67) 



which is negligible in the considered region. In particular, for the Gaussian beam the 
replacement (|67|) means the following: 



iV^exp [-{ul/cf] < Np 
(cf. dnni), and (jini))- This estimate leads to a large negative value of 



(6^ 



(69) 



and to 

{^2)± = -Np\{ma)±\\ (70) 

Just expression (fTUj) is a new "subtraction term" derived in [T7| ITHj. 

The mistake of ^] consists in replacement (jU7|) . This replacement is true for the 
region of small frequencies uol/c <^ 1 when J = Np and Z = —Np + J = Np{Np — 1) A^^ , 
but such a replacing is completely wrong in the considered region of large frequencies 
(|29|) . To show this, we perform numerical calculation of the sum J. For a given collision 
of beams we can consider a set of the longitudinal proton coordinates {za} as a set of 
random quantities with some distribution w{z). We assume below that 



w{z) 



(71) 
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Now the sum 



with 



0=1 



C-iS 



C = cos{qZa) , S = ^ sm{qza) , q = uj/c 



a=l 



is also the random quantity as well as 



J = C' + S' 



(72) 



(73) 



(74) 



The quantities C and 5" are the sums of large numbers of random items. Therefore, one 
can expect that they distribute in accordance with the normal law: 



dW 



1 



exp 



dC ^/2^F7^Ac 
where Ac = \/ (? — 1? and 



{C - iVpc) 
2A/p(Ac)2 



dW 

Is 



exp 



27r A/ps2 



^2 



2A/ps2 



(75) 



cos (g^) 



+00 



wi^z) cos{qz) dz = e ^"^ , s = sin (qz) = - 



(76) 



cos^ (qz) = I (^1 + 



= sin^ (qz) = | ^1 — e~ 



Numerical calculations had been performed using the generator of random numbers from 
MATLAB. These calculations confirm the above distributions fl75p . In particular, it can be 
seen from Fig. 4 and Fig. 5 where the results of numerical calculations for ql = 1, Np = 10^ 
and ql = 10, Np = 10^ are presented for lO'' various sets of {za}- 
Since 

{C) = NpC, {{C-Npcf)=Np(^-c^), (^) = 0, {S'')=Np'^, (77) 

we have 

(J) = Np + Np {Np -l)c\ {Z) = Np [Np - 1) c2 ^ Nl exp [-{ujl/cf] . (78) 

Moreover, taking into account that in the considered case (Npc)^ <^ Np, we find that (in 
contrast to (jHH|) . ^ ) 

{J) = Np, \{Z)\<^Np. (79) 
The distribution of the random quantity J becomes very simple at gZ ^ 1: 



dW _ 1 
d7 ~ iV^ 



-J/Np 



^0) 



The results of numerical calculations, presented on Fig. 6 , confirm Eq. ()8()|1 . It should 
be noted that distribution (jHIH) is rather wide. 



A J = ^{J^) - (J)2 = N, 



p ' 



(81) 
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therefore, the averaging over many various sets of {za} is necessary to obtain the stable 
result for (J). 

This consideration shows that the effect, derived in [T71 HH], is absent just in the 
region, discussed in these papers. At the end of this section we reconsider the experiments 
analyzed in paper [T7j . 

The HERA experiment jT2]- In this case Ee = 27.5 GeV and / = 8.5 cm, therefore, 
Ec = 27 keV. For the observed photon energies E^ = 2 ^ 8 GeV, the parameter 

--^>10\ (82) 

and the coherent contribution is completely negligible. Therefore, the new correction to 
the previous results on the level of 10 %, obtained in |17^, is, in fact, absent. 

The VEPP-4 experiment PJ. In this case Eg = 1.84 GeV and I = 3 cm, therefore, 
Ec = 0.34 keV. For the observed photon energies E^ > 1 MeV, the parameter 

^r.^>10^ (83) 

and the coherent contribution is completely negligible. 

The case of a "typical linear collider" with Ee = 500 GeV and E^ = E'e/lOOO. This 
example, considered in paper ^7j, is irrelevant for the discussed problem, since the co- 
herent radiation (called in this case heamstmhlung) at a typical linear collider absolutely 
dominates in this very region over the ordinary incoherent bremsstrahlung — see, for 
example, the TESLA project [22] and Sect. 3 of 



7 CONCLUDING REMARKS 



In the present paper we had performed analysis of the coherent and incoherent contribu- 
tions to the bremsstrahlung in conditions ()29j) when the incoherent contribution dominates 
but large impact parameters give an essential contribution to the standard cross section. 
In this conditions the known correction ()51|) to the standard cross section is determined 
by the transverse distribution of particles in the beams. 

We take into account correlations of particles in the beam. The corresponding correc- 
tion to the standard cross section is given by Eq. ()52|) and it is determined by correlations 
of particles in the transverse as well as in longitudinal coordinates. 

Through the paper we consider MD-effect in bremsstrahlung. The MD-effect for the 
e^e~ pair production (for example, in the reaction 76 — ^ e^e~e) can be considered in the 
same manner — for detail see Sect 7.1 from review [2]. 

We had shown that replacing the sum J (j66|l by the integral (j67jl in conditions (j29|) is 
incorrect. As a consequence, papers [dllH] are incorrect as well. 
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APPENDIX. AVERAGE PRODUCT OF PARTICLE DENSITIES IN THE BEAM 



In the classical limit ()24|) the average product of particle densities in the proton beam is 
given by the double sum over protons in the beam 

S{T,r')= 5^(5(r-rJ5(r'-r,)). (84) 

a, 6=1 

We split this expression into sum with a = b and sum with a ^ b: 
Sir,T') = S^ + S2, ^i = (5(r-r') 5^(5(r-r„)), = ^((5(r - rj 5(r' - r,)) . (85) 

To perform the averaging, we introduce the average proton distribution function 

/(r) = n,{v)/N, (86) 

with the normalization 

j f{r) d^r = l. (87) 

It gives 

(5(r - r„)) = y 5(r - r„) /(r,) = /(r) . (88) 

Note, that quantity (5(r — r^)) does not depend on index a and, therefore. 

Si = 5(r - r') Np /(r) = 5(r - r') n,(r) . (89) 

If we can neglect correlations between the particle coordinates, then the average prod- 
uct (5(r — Fa) (5(r' — r?,)) for a 7^ 6 can be presented as the product of two averaged 
factors: 

(5(r - r,) 5(r' - r,)) = (5(r - r,)) (5(r' - r^)) for a ^ 6 . (90) 
As a consequence, 

S, = 5^(5(r - rj) {6{r' - r,)) = N,{N, - 1) /(r) /(r') . (91) 

If we do not neglect the correlations between the particle positions, we should introduce 
the correlation function C (r, r') as follows 

S2 = J2{S{r - r,)) {6{r' - r,)) + C(r, r') . (92) 

a^b 

In that case we obtain instead of (PT|) the expression 

^2 = N,iN, - 1) /(r) /(rO + C(r, r') . (93) 

As a result. 

Sir, r') = N,{N, - 1) /(r) /(r') + 5(r - r') ^^(r) + C(r, r') . (94) 

Since in right-hand-side of this equation the first and the second items usually do not 
compensate each other, we can use approximation 

N,iN, - 1) ^ iV^ (95) 

and, therefore, 

5(r, r') = ripir) np{r') + (5(r - r') ^^(r) + C(r, r') . (96) 
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Figure 4: The curves s and c are given for dW/ dS and dW/ dC, respectively, in accordance 
with Eqs. (|75p and (j76p for ql = 1, Np = 10^. The histograms represent results of 
numerical calculations for distribution of S and C defined in (f73|) for 10"^ various sets of 
random numbers {za/l} 
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Figure 5: The same as in Fig. 4, but for ql = 10, Np = 10^ (in this case the curves s and 
c are practically coincide). 
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Figure 6: The curve is given for dW/dJ in accordance with Eq. ()8()|1 for ql = 10, Np = 10^. 
The histogram represents the result of numerical calculations for distribution of J defined 
in ()66|) . (f7^ for 10^ various sets of random numbers {za/l}. 
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